Comment on "Phase ordering in chaotic map lattices with conserved dynamics" 
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Angelini, Pellicoro, and Stramaglia [Phys. Rev. E 60, R5021 (1999)] (APS) claim that the phase 
ordering of two-dimensional systems of sequentially-updated chaotic maps with conserved "order 
parameter" does not belong, for large regions of parameter space, to the expected universality class. 
We show here that these results are due to a slow crossover and that a careful treatment of the data 
yields normal dynamical scaling. Moreover, we construct better models, i.e. synchronously-updated 
coupled map lattices, which are exempt from these crossover effects, and allow for the first precise 
estimates of persistence exponents in this case. 



I. INTRODUCTION 



The dynamical phenomenon of domain growth occurs 
in many different physical contexts. Once a fairly well 
established subject Q, it has been recently the center 
of renewed interest for several distinct reasons: one is 
the advent of new quantifiers of the associated dynam- 
ical scaling regimes, such as first-passage or persistence 
exponents, i.e. the rate of algebraic decay of the proba- 
bility for a given point in space to have remained in the 
same phase since some initial time . Another reason is 
the natural question of the extent of known "universality 
classes" to new types of systems, e.g. spatio-temporally 
chaotic ones. 

In this context, the recent study of the ordering prop- 
erties of chaotic coupled map lattices (CMLs) possess- 
ing several symmetric phases in competition brought up 
some intriguing results Q : in the simple case of two com- 
peting phases and a non-conserved order parameter, the 
"normal" growth law L{t) ~ t^^^ with z — 2, where 
L{t) is the single lengthscale characterizing the coarsen- 
ing pattern, was observed but with some exponent z > 2 
continuously varying with parameters. However, this 
was later shown |^ to be only a (slow) transient behav- 
ior due to the non-trivial effect of space-discretization 
in these deterministic systems. For larger lattices and 
longer times than those considered in H, L was shown 
to behave normally when plotted against t^^^. 

In a recent Rapid Communication, Angelini, Pellicoro, 
and Stramagha (APS) |5j, motivated by the above study 
of non-conserved order parameter CMLs, presented a 
class of sequentially-updated lattices of chaotic maps de- 
signed to investigate the case where the order parame- 
ter is locally conserved. In this case, L is also expected 
to grow algebraically with time, but with z = 3 
APS claim, however, that larger exponents are commonly 
found. Here we show that APS were misled by their 
treatment of data and that in fact the normal (z = 3) 
growth law is observed in all cases. We argue, moreover, 
that fully deterministic, synchronously-updated, coupled 
map lattices which conserve the order parameter can 
be easily constructed following the ideas of Oono and 



Purl 01, and we show that these systems behave very 
smoothly, enabling the precise measurement of persis- 
tence exponents in this context. 



II. REVISITING APS RESULTS 
A. The hybrid map lattices of APS 

The lattices of maps introduced by APS are hybrid in 
several ways: a given local map / is first applied to all 
sites Xi of the lattice (a deterministic and synchronous 
operation), then pairs of nearest neighbors are sequen- 
tially and regularly visited and swapped probabilistically. 
(The regularity of the sweeps of the lattice is at the ori- 
gin of the anisotropy of the domains in Fig. 1 of [^.) 
Furthermore, such systems are not "coupled" maps as in 
usual CMLs, since the values taken by the sites are not 
influenced by the swaps (they are always taken according 
to the invariant measure of the local map). 

These systems are designed to mimick Ising systems 
(with fluctuating couplings) corresponding to the "spins" 
CTi = sg'n{xi). The local map / has not much impor- 
tance, and it is convenient to choose, following [^,D, 
an odd map of the [—1, 1] interval with two symmet- 
ric attractors.. The energy of one configuration is given 
hy E — —^(^ij-^XiXj where the sum is over nearest- 
neighbor pairs. The exchange probability reads Pgwap = 
1/(1 + exp(/3A£^) where (3 is the inverse temperature 
and Aii' is the energy change of the swap. The zero- 
temperature limit is deterministic: swaps are effective if 
and only if they decrease the energy. 



B. Domain growth is normal 

We have performed numerical simulations of the APS 
system at zero temperature with the piecewise linear lo- 
cal map used in H: 



liX if X e [-1/3,1/3] 

f{x) = { 2Ai/3 -^lX if a: G [1/3, 1] 

-2^i/'d- ^ix if a: e [-1,-1/3] 



(1) 
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with fi = 1.9. Coarsening occurs, with, again, a strong 
anisotropy due to the mode of update. The growth of 
L, defined as the width at mid-height of the two-point 
autocorrelation function, is slow at short times, but then 
reaches the expected t^^^ behavior (Fig. |l]a), contrary 
to the claims of APS. The short-time behavior may be 
mistaken for anomalously slow algebraic growth (with an 
exponent close to the value 1/z = 0.07 reported by APS) 
when logarithmic scales are used, but a closer inspec- 
tion shows a systematic increase of the local exponent 
(Fig. |l|b). As a matter of fact, the system is so anisotropic 
that domains elongate in time (Fig. |l|a). This is due to 
our choice of updating x-wise pairs before y-wise pairs. 
Alternating this order would presumably suppress this 
effect. 
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FIG. 1. Domain growth in the APS system with the local 
map (|l|) at zero temperature (/3 = oo). Lattice of size 1024^ 
with periodic boundary conditions. Solid lines: L defined as 
the width at mid-height of the two-point correlation function 
estimated along the a;-axis of the lattice. Dashed lines: same 
but along the y-axis of the lattice, (a) L vs t^^^\ insert: local 
slopes, (b) ln(L) vs ln(t); insert: local exponent. 

Runs of the same system at finite temperatures indi- 
cate that domain growth is faster and that the crossover 
to the z — behavior occurs sooner. We are confident 
that similar results hold for the complicated local map 
mostly used by APS. As a conclusion. Fig. 3 of |^ has 
to be replaced by the variation of the prefactor of the 
L oc t^^^ law, similarly to the final conclusions of [Q for 
the non-conserved order-parameter case. 



C. Persistence scaling is hard to measure 

The persistence probability p{t) = Prob{(Ti(t') — 
(^i{to), Vi' G [^0,^]} is usually observed to decay alge- 
braically with time {p (t/to)^^) in systems with alge- 
braic growth laws. But persistence scaling for conserved 
order parameter systems is notoriously difficult to ob- 
serve Q. An additional difficulty lies in the fact that 
the available models only show coarsening at finite tem- 
peratures, so that one has to resort to block-scaling of 
persistence. This is also the case of APS systems, even 
at zero temperature, since the chaotic fiuctuations of the 
"couplings" amount to a finite temperature. This, by the 



way, is the reason why APS systems coarsen in this case. 

Needless to say, the estimates of persistence exponents 
presented in |^ are then highly unreliable, if only because 
of the slow crossover for the growth law of L. Ideally, 
since persistence is a complex quantity involving all times 
since the reference time to , one should in principle choose 
to in the asymptotic scaling regime and simulate the sys- 
tem up to times t ^ tg- Given the typical values of the 
crossover times (Fig. 0) this is hardly possible. Another 
difficulty for the APS systems is the possible infiuence 
of their strong anisotropy on the persistence exponent 9 

Rather than trying to measure properly persistence 
scaling in APS systems, a possible but difficult task, we 
now turn ourselves to truly deterministic models, i.e. reg- 
ular coupled map lattices, which are devoid of the draw- 
backs underlined above for APS systems. 



III. WELL-BEHAVED DETERMINISTIC 
MODELS 

A. Oono-Puri style CMLs 

Deterministic models for phase ordering of conserved 
systems were introduced by Oono and Puri Q . Dynam- 
ical scaling was observed to hold for these CMLs, with 
a growth law compatible with the expected z = 3. We 
now present similar models which, in addition, can be 
constructed for any local map. 

A usual CML, such as that studied in |^,^ for the non- 
conserved order parameter case, can be written: 



(2) 



where g is the coupling strength, TV is the number of 
neighbors in the choosen coupling range, and the sum is 
over these neighbors. 

Following H, a CML conserving exactly the continu- 
ous local field can easily be constructed as: 



-y 



(3) 



where the last term corresponds to the extra Laplacian 
in the Cahn-Hilliard equation. Conservation of the con- 
tinuous field X is obvious. On the other hand, strictly 
speaking, the discrete field a — sgn(x) is not exactly 
conserved and fluctuates slightly, because the last term 
in ^ may change the sign of sites situated in domain 
walls. 

The synchronous mode of update prevents excessive 
(i.e. other than lattice-derivated) anisotropy. The above 
structure insures that "true" zero-temperature regimes 
are observed if the local map possesses two disjoint at- 
tractors. Changing the nature of these attractors (fixed 
points, limit cycles, chaotic sets), one can study compe- 
tition between phases of different nature. 
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B. An example 



IV. CONCLUSION 



We now present results obtained on a particular case 
of the models defined above. More comprehensive results 
will be reported elsewhere For simplicity reasons, 

we again choose the map given by Eq. ([T]). For "extra 
smoothness" , the Moore neighborhood on the square lat- 
tice (A/" = 8 neighbors of equal weight) was used. 

As in domain walls are strictly pinned for 

small g. For too strong coupling, on the other hand, 
antiferromagnetic-like phases appear. There is, how- 
ever, an intermediate range of g values for which domain 
growth proceeds forever between two weakly-chaotic 
phases. The expected 2 = 3 law is then easily observed 
even at short times and in log- log scales (Fig. 0a). 
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FIG. 2. Phase ordering and persistence scaling in a 
Oono-Puri style CML. Lattice of size 8192^ with periodic 
boundary conditions, (a): log(L) vs log(t); insert: local ex- 
ponent (The slightly smaller values of l/z recorded can be 
shown to be due to space discretisation effects.) (b): log(p) 
vs log(t) for various initial times to. From bottom to top: 
to = 0, 10, 10^, and 10^. 

The above CMLs reveal their strongest advantage 
when persistence scaling is considered. As already no- 
ticed, they show normal coarsening with true zero tem- 
perature. This allows to avoid studying block persistence 
scaling, a somewhat tedious task. Figures shows per- 
sistence decay for different reference times to. Nice scal- 
ing is easily observed. This constitutes, to our knowledge, 
the first clean evidence of algebraic decay of persistence 
in a two-dimensional conserved order parameter system. 
Our results give an exponent 6 ~ 0.25(2), i.e. a value 
larger than that observed for the non-conserved order 
parameter case (for which 9 ~ 0.20 - 0.22 ||,|l|l). 

The above CMLs constitute an excellent base for a re- 
liable study of persistence scaling in conserved order pa- 
rameter systems. Ongoing work is probing the degree of 
universality of the persistence exponent measured above, 
alongside a similar study for the non-conserved order pa- 
rameter case for which this issue is still unresolved MfiM ■ 



In this Comment, we showed how APS were misled 
in their interpretation of simulation data and that their 
conclusions about possible non-trivial values of the dy- 
namical exponent z in chaotic systems with conserved 
order parameter dynamics do not hold. To a large ex- 
tent, these systems can be seen as too close to the Ising 
model with Kawasaki dynamics, which is well-known to 
be difficult to study numerically (although z = 3 scaling 
is now well documented [p^). 

We introduced a class of CMLs which are devoid of 
all the problems encountered in APS systems and which 
show normal scaling already at early times. Moreover, 
these systems also present nice scaling behavior for the 
persistence probability, whereas similar investigations in 
APS systems are riddled with problems. We showed 
unambiguously that two-dimensional systems with con- 
served order parameter domain growth show an algebraic 
decay of persistence. Future work will try to assess the 
universality of both the so-called Fisher-Huse exponent 
and the persistence exponent 9 in such systems. 
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